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1. Introduction 

The Gardner equation appeared a long time ago when Miura [1, 2] introduced the 
remarkable transformation 

u = v"^ + v^, (1) 

connecting the Kortewcg-dc Vries (KdV) equation, 4ut = u^x — Quu^, to the modified 
KdV (mKdV) equation, Avt = v^x — Qv'^Vx- Both equations are ideal prototypes for 
integrablc models. The transformation (1) is highly nontrivial and relates solutions of 
the two nonlinear equations. 

Miura [1] pointed out that the KdV equation is invariant under Galilean 
transformation. In fact, if we change the coordinates according to 

x^^ + ar, t ^ T, a = const., (2) 

then ^ = ^5 §i ~ §^ ~ ^§1 ^"^^ KdV equation becomes Aur — Aau^ = — 6uu^. 
The undesirable —Aau^ term on the left hand side is canceled by the nonlinear term 
—6uUx on the right hand side, by a constant shift in the field variable. Therefore, the 
KdV equation is invariant under the transformation (2) with 

u{x,t) ^ u{^,t) + '^a. (3) 

Following the same idea for the mKdV equation we get Avr — Aav^ = v-^^ — Qv^v^. Because 
of the quadratic term in the nonlincarity. a constant shift in the field can still remove the 
undesirable —Aav^ term, but at the cost of an extra nonlinear one. Then, the mKdV 
equation is not invariant under the transformation (2) but it is transformed into the 
Gardner equation 

Avr = — 6^^^^ — 12jj,vv^, jj, = , (4) 

if we shift the field 

v{x,t)^v{^,T) + {lay^\ (5) 

Equation (4) is evidently interesting because it may be considered as combination of 
the KdV and mKdV equations or as an integrable deformation of the KdV or mKdV 
equations. In a sequence of three papers, Wadati [3, 4, 5] obtained the n-soliton solution 
using the inverse scattering transform, Backlund transformation and showed that (4) is 
indeed a completely integrable Hamiltonian system. 

It is known that integrable soliton equations appear as members of a more general 
structure called integrable hierarchies. In [6] the results previously obtained for the KdV 
equation individually were generalized to the whole KdV hierarchy and its deformation 
related to the Gardner transformation was discussed. 

The properties of the Gardner equation is fairly well known and can be obtained 
through the properties of the mKdV equation, and is subject of current research [7, 8]. 
Besides the mathematical interest, equation (4) has a wide application in atmospheric 
and ocean waves and was extensively studied by Grimshaw et al [9] where they obtained 
the large-amplitude table-top solitons, often observed in ocean coastal zones, and 
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breathers [10]. Recently, solutions of the Gardner equation were used to construct static 
solutions of the Gross-Pitaevskii equation, known to describe dynamics of Bose-Einstein 
condensates [11]. 

A general algebraic approach for the construction of integrable hierarchies and its 
solitonic solutions can be formulated in terms of representation of affine Lie algebras. 
The sohton solutions can be systematically constructed by the dressing method, which 
connects a trivial (vanishing) vacuum to a nontrivial soliton solution employing vertex 
operators. See for instance [12, 13, 14, 15, 16] and references therein. 

The transformation (5), connecting solutions of the Gardner and mKdV equations, 
clearly shows that both solutions have different boundary values due to a constant shift 
in the field variable. This implies that both equations have different vacuum solutions. 
In particular, those solutions of the Gardner hierarchy with trivial vacuum are related 
to nonvanishing vacuum solutions of mKdV equation and vice-versa. 

In [17] we have extended the mKdV hierarchy for negative flows by introducing 
negative even graded Lax pairs. The nonlinear equations within that class do not 
present solutions with vanishing boundary condition and, in order to circumvent this 
problem, a deformed vertex operator was introduced. The present paper is a natural 
continuation of [17] in which the same algebraic formulation is employed to construct 
solutions of the Gardner equation. Interesting solutions arises where a table-top soliton 
or a kink can coexist with normal solitons. We also obtain dark solitons, breathers and 
the wobble [18, 19] solutions. 

In sections 2 and 3 we briefly introduce the KdV and mKdV hierarchies, 
respectively, and show that Miura transformation represents a map between both 
hierarchies and not only between two single equations. Our main results are contained in 
section 4 where we propose a deformed hierarchy grounded on si2 Kac-Moody algebra, 
containing the Gardner equation as one of its members. In section 5 we introduce a new 
vertex operator and in section 6, construct explicit solutions for the whole hierarchy. 
Goncluding remarks are in section 7. 

2. The KdV hierarchy 

Following [20] , the mathematical construction of the KdV hierarchy can be introduced 
by the Lax equation in terms of pseudo-differential operators 



where L = — u{x, t) is the Schrodinger operator. The subscript L_|_ denotes the 
differential part of the operator. We then calculate 



Lt = [Lf,L] 



n = 1, 3, 5, . . . 



(6) 



U2x) d - ^U3^ + fuu, 



(7) 
(8) 



Substituting (7) in (6) we obtain the KdV equation 



6uu. 



(9) 
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Repeating the same calculation with (8), we have the Sawada-Kotera equation^ 

16Ut — U^x — lOuUsx — '^0UxU2x + SOu^Ux- (10) 

Higher order equations are obtained from (6) and its soliton solutions are well known 
[20]. It is important to note that while (9) is invariant under Galilean transformation, 
(10) and higher order KdV equations are not. 



3. The mKdV hierarchy 



The mKdV hierarchy can be constructed from a zero curvature condition based on the 
affine Kac-Moody algebra 3^2 [12, 13, 14, 15, 17], generated by 

si2 = {E^:^\ (11) 

The s£2 



together with the spectral derivative operator d and the central term c 
commutation relations are 



77, (n) 77. (m) 



= 0, 



where T^") G s£o. 



[12) 



G s^2- The operator Q = \H^^^ + 2d defines the grading operation 
[(5,T'^")] = nT^"\ that decomposes the algebra into even and odd graded subspaces 



si 



(2n) 



(2n+l) 



(13) 



Considering this grading structure, it follows from the Jacobi identity that for T^*) G si. 
and T(J) G sY^^ we have [T«,T(^)] G s/^^^\ 



(0 



The mKdV hierarchy is then defined by the zero curvature equation 

dx + Ej^^ + E^_^l + vH^^\ dt + D^""^ + D^^-^) + • • • + D^^^ 







(14) 



where n — 1,3,5, . . ., v — v(x,t) and D^^^ G I^ the construction of integrable 

models from the zero curvature equation, only the loop-algebra, that corresponds to set 
c = in the commutation relations (12) is employed. Equation (14) is solved grade 
by grade, starting from the highest one, until the zero grade projection leads to the 
nonlinear time evolution equation. This procedure works in the following way. Let 
n = 3 for example, then from the grading structure (13) the operators involved in the 
construction (14) must be linear combination of the algebra generators in the form 

D^'^ = a,E^^ + hE^^\ = C2H^'\ (15) 

= aii?(°) + h-,E^^\ = coH^'\ (16) 

The coefficients {ai,bi,Ci) will be determined as functions of the field v{x,t), by 
projecting the zero curvature equation (14) into each graded subspace. In this way, 
wc obtain the nonlinear partial differential equation, from the zero grade projection, 

X The coefficients of this 5-th order KdV are different from the original Sawada-Kotera but the nonhnear 
terms are the same. 
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and also its Lax pair. Note that only the commutation relations (12) (with c = 0) are 
used in this calculation and no matrix representation is needed. So, for n = 3 we get 
the well known mKdV equation 

4:Vt = - dv'^v^ (17) 

while for n = 5 we find 

IQvt = V,, - lOv^sx - ^Ovv,V2. - lOvl + 30v^v,. (18) 

Without loss of generality in the above calculations, all integration constants were chosen 
to vanish. 

Using the Miura transformation (1) in (9) and (10), we can relate them with 
the corresponding equations in the mKdV hierarchy, (17) and (18), respectively. This 
operation can be expressed for both equations as 

KdV {u^v^ + v^) = {2v + d^) mKdV (v) . (19) 

After verifying this expression for this two particular cases, wc now show that, in fact, 
the Miura transformation holds for all higher orders, i.e. the Miura transformation 
is a map between solutions of the entire mKdV hierarchy into solutions of the KdV 
hierarchy. Consecutive time evolutions of the KdV hierarchy are given by the recursion 
formula 

«W2=-R^t„> Ut,=u^, (20) 

where R = ^d^ — u — \uxd~^. Substituting u = v"^ + Vx both sides of (20) we can 
write it in a factorized form as 

«t„+2 = (2^ + dx) vt„^, = {2v + dx) {{dl -v"- Vxd-^v) v^^. (21) 
Define now R! = \d'l — — Vxd~^v, thus 

Vt„+2^R'Vtn, Vt.^Vx, (22) 

which shows that entire mKdV hierarchy is recursively generated by R' and (19) holds 
for n — 1,3,5, 



4. A hiercirchy containing the Gardner equation 

Motivated by the transformation (5), we now propose a deformation of the mKdV 
hierarchy (14), 

dx + E^^^ + E^H + (Ai + v)H^''\dt + + D^^-^^ + • • • + L>(°)] = (23) 

where // = const., v — v{x,t) and n — 1,3,5, In the usual algebraic construction 

[14], as used in (14), the semi-simple element E = Ea^ + E^}\ is responsible for 
the algebra decomposition s^2 = IC Q) Ai, its kernel JC and image Ai subspaces. In 
(23) the deformed element Ea^ + E'^1, + ^H^^^ contains both /C and M. components. 
Obviously, the modification from (14) to (23) corresponds to a simple translation in the 
field V ^ V + 11, however, it introduces important changes in constructing its solutions 
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through the dressing method. Reconsidering (14), but now with a nonvanishing constant 
boundary condition, v Vq, the situation is quite the same as in (23) with f — )■ 0. 
Therefore, solutions with vanishing boundary condition of (23) are related to solutions 
with nonvanishing boundary condition of (14). 

Solving (23) for n = 3, using the same procedure described in the previous section 
for the mKdV hierarchy, we get 

vt = -^v^x - [l^ ~ 2 / ~ ^l^vv^ - -V (24) 

where a is an arbitrary integration constant. Choosing a — 2//^ we have the well known 
Gardner equation 

^Vt = vzx - i'^IJ'VV^ - 6v\^ (25) 
and its Lax pair§ 

A = + E^_}1 + {fi + v) (26a) 
B = + + {f, + v) H^'^ + l{v,-v^- 2f,v + 2f,') - 

- i (^^0. + + 2iiv - 2//) E^ll + \ {v2x - - 6iiv^ + 4//) H^""^ . (266) 
Solving now (23) for n = 5 we obtain 

16Vt = fsa; - lOv^s^ - A0vv^V2x - lOf^ + SOv%^ - 

- 20lIVV3a: " 40//Va,V2a; + 12^11^ v'^V^ + 120//^^. (27) 

The two arbitrary integration constants that show up are conveniently chosen: a = 4/i^, 
/3 — 6/x^. Equation (27) is a combination of (18) and (10). Note that when // — >■ 0, (27) 
becomes (18) in the same way that (25) becomes (17). We point out that (27) is not 
obtained from (18) by the Gahlean transformation (5). The Lax pair for (27) is 

A = E^^^ + E^ll + {ii + v) H^'\ (28a) 
B = + E^^l + ifi + v) + 1 (^^ _ ^2 _ ^ 4^2^) ^(1) _ 

-l{v^ + v^ + 2i2v - Ajj^) E^^l + \ (t;2, - 2v'^ - Q^v^ + A^^v + 8/^2) H^^^ + 

+ \ {Vzx - Gv'^V,, - 2vV2x +vl + 3v'^ - UflVVx - 2lJV2x + 

+Ai/vx + I2fiv^ + S^i^v^ - Sfi^v + 8fi^) - 

- I {V3x - Qv^Vx + 2vV2x - vl - 3v'^ - UflVVx + 2/i?;2x + 

+4/i\^ - 12/xw^ - + 8iJ,\ - 8/i^) e'lI + 

+ {v4x - Wv\2x - Wvvl + - 20fxvv2x - I0f^vl+ 

+30//^;^ + AOiJ^v^ + 16/x^) H^'^^ . (286) 

The hierarchy defined in (23) can also be recursively generated through the following 
pseudo-differential operator 

Vtr.+2^RVt„, Vt^^Vx, (29) 

R = Idl -{v + 2^^)v- V, {ixd-' + d-H) . (30) 
§ L = d + A,Lt = [L,B]. 
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5. Vertex operator 

The usual $£2 vertex operator with principal gradation 

00 

v{k)^ j2 [^^"^ - ¥noc + ^-'E<^:) - K-^E'f!:'^] , (31) 

n=— 00 

solves the mKdV and sinh-Gordon equations with vanishing boundary condition. Here 
we introduce a modified vertex operator that generalizes and takes into account the 
nonvanishing boundary value problem for the mKdV hierarchy and henceforth, the 
Gardner hierarchy with vanishing boundary condition. 

Define the following vertex operator depending on the parameters {ni, /J^) 



2\ 



(32) 



This vertex was proposed in [17] to solve the negative even grade part of the mKdV 
hierarchy. Note that when — > the vertex (31) is recovered. The parameter fi, also 
present in equation (25), is related to the nonvanishing boundary condition of the mKdV 
hierarchy through the identification fj, ■v^ vq, where v —> Vq is the constant value of the 
field in \x\ 00. Consider the operator 

Qm = Ej^^ + Et+'^ + (33) 

which corresponds to the vacuum v = configuration of the Lax component A^ac = ^0 
in (26a), for m = 0. Taking the commutator of (33) with the vertex (32) we verify that 

[flm, Vi\ = -2Ki {k^ - ■ (34) 

As will be clear in the next section, (34) determines the dispersion relation for all 
nonlinear evolution equations in the hierarchy. Consider the highest weight states for 
s£2 — {|Ao), |Ai)} — which obey the following actions: Ea^\\j) = 0, E^^\Xj) = and 
ii-W|A/) = forn > 0, //(°^|Aj) = 5ji\\j) and c\\j) = for j = 0, 1. The adjoint 

relations are: — H^~"'\ E^^ — E^~^^ and — c. Taking (32) between the highest 

weight states, 

{\j\Vi\Xj) = ^i±^ where (Tq = -1, ai = 1. (35) 
A more involved calculation of two vertices yields 

{Xj\ViVi\Xj) = ^ '—^ -^CLij where Qij = — - — . (36) 

4:KiKl \Ki -\- Ki J 

It is possible to prove that the expectation value of a product of n vertices decomposes 
into 

n n n 

i=l i=l i,k=l,i<k 

The matrix elements (36) and (37) determine the nonlinear interaction between solitons. 
Notice that the nilpotency property of the vertex representation is a direct consequence 
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of (36) when Ki = ki. This corresponds to the physical interpretation that when 
Ki — )■ Ki two sohtons degenerate into a single soliton (no interaction). Equation (37) 
also corresponds the physical property that n solitons interact in pairs. 



6. Dressing the vacuum 

The two dressing group elements 0± correspond to gauge transformations mapping 
trivial vacuum potentials, Aq and Bq, into nontrivial ones involving field dependent 
potentials, A and B [15, 13] 

A = e±Ao©i^ - d^e±ei\ (38a) 

B = e±Soe±^ - dte±ei\ (38&) 

They are factorized into positive and negative grade generators 

0+ = e^'^^^^.. e_ = e^'-\^'-^ . . (39) 

which, together with the zero curvature condition yields, 

0:^0+ = ^0^*0^ (40) 

where G is an arbitrary constant group element and ^'o is such that Aq — — ^ 
and Bq — —dt^o'^Q^. For the Gardner equation (25), taking in (26a) and (26&) we end 
up with the vacuum potentials Aq — and Bq — + /i^^q, therefore 

*o = exp [-Qox - (Qi + /i^l^o) t] (41) 

where flm is defined in (33). The vacuum takes into account information about the 
boundary condition at |a;| — )■ oo (t; — )■ 0). 

The dressing method requires the use of highest weight states representation so it 
is necessary to include the central term c in the calculations. In the construction of the 
models, because c commutes with every other operator, the zero curvature equation is 
invariant under addition of a central term. Therefore, we change A ^ A ~ u^c, where 

is one of the Lax operators and z/ is a function that will be 
determined. Solving for the zero grade projection of (38a) with 6+ we find 

exp [X^^^] = exp [^ii-^o) + i/c] , -u = -cf)^. (42) 

Using this result when taking the left hand side of (40) between highest weight states, 

e'^ = (Ao|e:ie+|Ao) = tq, e<^+^ = (Aiiei^e+lAi) = n, v = d^ln^. (43) 

Choosing now G — Y[i=i projecting the right hand side of (40) between these 

states 



(A,inexp{^oV^^^o'}|A, 



1=1 
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= (A,inexp{e«»(-'V,}|A,) (44a) 

i=l 

n 

1=1 

where in (44a) we have used the fact that Vi — V{K,i) are eigenstates of with 
eigenvalues given by (34). It therefore follows the dispersion relation, 

= 2KiX + {2Ki [k] - /i^) + 2Kiix^] t = 2KiX + 2K^^t (45) 

where in (446) the nilpotency property of the vertex Vi was enforced. The general 
formula (446) is valid for all equations within the hierarchy and can be solved explicitly 
by using the factorization property of the vertices (37) and (36). The explicit space-time 
dependence, however, is specified according to the choice of vacuum potentials Aq and 
Bq for each individual model. Let us consider for instance the vacuum for potentials 
(28a) and (286), Aq = VLq and ^0 = ^2 + 2ii^Vti + ii^VLq. From (34) we therefore have 

= 2KiX + 2Klt. (46) 
6.1. Nonvanishing boundary condition for the mKdV hierarchy 

We have constructed solutions to the hierarchy associated with the Gardner equation 
where — >■ when |a;| 00. We now address the nonvaninshing boundary value 
problem for the mKdV hierarchy. 

Recall that in (23) if we take = we recover the mKdV hierarchy (14). When 
taking a constant vacuum v ^ vq m. (14), the term vq plays the role of and hence 
throughout this subsection you shall consider all expressions after (32) with // replaced 
by Vq. 

Reconsider the zero grade projection of (38a) for 0+, but now with Ao = + 
E^ll + vqH^^^ and A = E^a^ + E^^l + v{x,t)H^^^ - u^c. The analogue of (42) becomes 

exp = exp [{voX + 0) H^°^ + uc] , -0^, (47) 

and (43) implies the following nonvanishing boundary solution 

V ^vo + d^\n-. (48) 

Equation (446) docs not change except for the dispersion relation. The vacuum 
potentials arc obtained by setting v Vq in (14). For the mKdV equation (17) we 
have ^0 = ^0, -Bo = ^1 — ^vo^o and for (18) Ao = VIq, Bq = Q.2 — \vl^i + ft'o^o- Then, 
using (34), we have the respective dispersion relations 

= 2KiX + {2kI - ?,vIk,) t, (49) 
= 2KiX + [24 - ^vIkI + f t;^fi:,) t. (50) 

Note that the dispersion relations depend on the boundary condition, so the velocity of 
this solitons depend on the vacuum field Vq. 
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(a) (b) (c) 

Figure 1: Soliton for the Gardner equation when ki > ^ and /i. > 0. We choose fx — 2 
and Ki — (3, 2.0001, 2) in figures (a, b, c), respectively. When ki — >■ the solution 
shown in (a) deforms, as shown in (b), until it becomes a kink when ki — /i as shown 
in (c). 



6.2. Examples 

We now present some explicit examples of the exact solutions we have obtained. The tau 
functions for a general n-soliton solution is given by (446). The solution is then obtained 
from (43) or (48), with space-time dependence specified by the dispersion relations (45)- 
(46) or (49)-(50), according to the Gardner or mKdV hierarchies, respectively. 
The general expression for 1-soliton solution is 

Tj = l + (\/i),e«^ (51) 

where j = 0, 1 and = {Xj\ • |Aj). The vertex clement is given by (35) with n ^ Vq 
for the mKdV hierarchy. The 2- and 3-soliton solutions are respectively given by 

r, = 1 + {V,),e^^ + {V,)je^^ + (K^),e«^+«^ (52) 

and 

rj^l + {V,)je^^ + {V2)je^^ + {Vs)je^^ + {V^V,) je^^+^- + 

+ (T4T^3),e«i+«^ + {V2Vs)je^'+^' + {ViV2)j{ViVs)j{V2V3)je^'+^'+^\ (53) 

The 1-sohton solution of the Gardner equation, with dispersion (45), has different 
behaviour according to the values of Ki and /i. When Ki > /i and /i > the solutions 
are plotted in figure 1. Note that in figure Ic we have a kink (anti-kink if you are 
used to sine-Gordon terminology). When ni < but \k,i\ < // we have figure 2. For 

> n we recover the behaviour of figure la. The more interesting solution occurs 
when < Ki < /I, having the usual KdV solitons and also a maximum amplitude table- 
top soliton, as shown in figure 3. The same graphs apply to (27), the only difi^erence 
is in the dispersion relation (46). If we change sign n — )■ —fi, Ki — )> —Hi the graphs 
are reflected through y and x axis, so we have an elevation instead of a depression 
wave, figure 3b. In figure 3c we have the mKdV nonvanishing boundary solutions with 
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-1.5 -1 -0.5 0.5 1 -3 -2 -1 l 




(a) (b) (c) 

Figure 2: Soliton for the Gardner equation when ki < and < /x. /x = 2 and 
Ki — (—1.5, —1.9999, —2) in figures (a, b, c), respectively. 




(a) (b) (c) 

Figure 3: Sohton for Gardner (a, b) and mKdV (c) when < ni < fx. The usual 
KdV like soliton increases its amplitude and becomes narrower as increases, but 
when Ki — > ii~ it becomes the maximum amplitude table-top soliton. /i — 2 and 
K,i — (1.5, 1.9, 1.999999999) in (a). In (b) we changed sign fi — > — k,i — > —ki with 
the same numerical values. In (c) wc have mKdV dark solitons and table-top soltion, 
with the same numerical values as in (a). Note the background corresponding to Vq — 2. 

dispersion (49). The mKdV equation can also have a kink, hke in figure la, but with 
the asymptotes in +2 and —2 corresponding to the constant background Vq — 2. The 
same solutions also apply to (18) with dispersion (50). This solutions with a constant 
background are known as dark solitons [21, 22]. Dark solitons of the NLS equation has 
wide application in nonlinear optics. 

Substituting the previous solutions of the mKdV hierarchy in the Miura 
transformation (1) we obtain dark solitons of the KdV hierarchy. Unhke the mKdV, 
the KdV hierarchy does not have kinks or table-top solitons, therefore, the interesting 
solutions of the Gardner equation is inherited from the nonvanishing boundary solution 
of the mKdV equation. 
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(a) (b) (c) 

Figure 4: (a) breather for the Gardner equation, (b) breather for the mKdV equation: 
ljL — 2,a — 0.5, /3 — 5. (c) breather for the KdV equation: fj, — 2, a — 1.5, /3 — A. 




(a) i = 5 {h)t = 2 (c) t = (d) t = -5 

Figure 5: The waves travel to the left so you should read this figure from (d) to (a). 
Wc arc plotting a 3-soliton solution of the Gardner equation with parameters /j, — —2, 
Ki = -1.99999, Av2 = -1.5, K3 = -1.2. 



For the 2-soliton solution (52) we can have a combination of two usual solitons, one 
soliton with a kink or one soliton with a table-top soliton. Moreover, we can have the 
interesting breather, a spatially localized but oscillating solution, by choosing complex 
conjugate wave numbers Ki = a + if3 and k,2 = a — i(5 (the explicit expression for 
the breather is given by (A. 4) in Appendix A). Figure 4 shows the breather for the 
Gardner, mKdV and KdV equations. The nonvanishing boundary breather of the KdV 
is obtained through Miura transformation. 

In figure 5 we have a specific situation of the 3-soliton solution (53) . Note that the 
final profile is different from the initial one. The waves did not keep their initial form 
after interaction. The mKdV equation can also have the same kind of 3-soliton solution 
with a background vq. 

In [18] Kalbermann proposed the wobble solution for the sine-Gordon equation then 
Ferreira et al [19] showed that the wobble is a 3-soliton solution, where two solitons 
combine to form a breather and the third one is a kink. We follow this line of thought 
and take our 3-soliton solution (53) with ki — a -\- i[3, K2 — a — ij3 and = to 
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(a) (b) 

Figure 6: In (a) we have the wobble sohition of the Gardner equation with parameters 
Ki = ^ + 6i, /t2 = ^ — 6i, = /X = 2, t = —0.02. In (b) we have wobble for the mKdV 
equation: ki = ^ + 5i, K2 = \ — 5i, k,^ = vo = 2, t = —0.02. 

obtain the wobble solution for the Gardner equation (the explicit expression is (A. 9) in 
Appendix A). We can also have the wobble for the niKdV equation with k,^ = vq. This 
solutions are shown in figure 6. Despite the KdV equation can have a breather, it can 
not have a wobble because it docs not have a kink. We could also combine the breather 
with a table-top soliton and with usual solitons as Grimshaw et al considered in [10]. 

7. Concluding renicirks 

We have constructed an intcgrable hierarchy (23) that contains the Gardner equation 
(25) as one of its members. This construction is based on the Kac-Moody algebra 
s^2 with principal gradation. Besides the Gardner equation, another fifth order PDE 
(27) that is a combination of Sawada-Kotera and fifth order mKdV (18) was explicitly 
considered. 

We have introduced a new vertex operator (32) that used within the dressing 
method, enabled us to obtain explicit n-soliton solutions of the Gardner and mKdV 
hierarchies, this last one with nonvanishing boundary condition. Besides usual KdV like 
solitons our solutions contemplate table-top solitons and kinks. The 2-soliton solution 
can be used to form a breather and using the Miura transformation we obtained a 
breather for the KdV equation, also with nonvanishing boundary condition. 

The 3-soliton solution showed that the interaction of individual waves when a 
kink or table-top soliton is present changes its initial profile such that energy is always 
conserved. Combining a breather with a kink we also obtained the wobble solution for 
the Gardner and mKdV hierarchies. 

We stress that the our solutions are valid for the whole hierarchy of nonlinear 
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equations, the only modification relying on the dispersion relations. 

Further exploration of our results can be made by considering the nonvanishing 
boundary value problem for the non-abelian AKNS hierarchy that contains the NLS 
equation which has practical applications in nonlinear optics and water waves. 
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Appendix A. Breather and wobble solutions 



The general breather expression is obtained from (52) by setting complex conjugate 
wave numbers, ki — K2 — + i^, which implies that the dispersion relations will be 
in the form — — V ''<' Ki ^r some real functions 77 = r]{x, t) and C, — Ci^i t) to be 
determined latter. Prom (35) we have {Vi)j — {V2)* — aj — ib where 



2 (a2 + /32; 
2(a2 + /32)' 



+ 



((7o = -1, (Ti = 1), 



and from (36) 



{V^V2),^-^^{a^ + b^). 
Therefore, the general tau functions for the breather are 



1 + 2e'' 



Qj cos( + bsm( — 



2o? 



ia^ + b^) 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



The wobble is obtained from (53) with Ki = ^2 = « + i/3 and k,^ = fi, so 
^1 = Q = f] ^ K ci'iid ^3 = fl^J. '^^ ^ r^s-l function depending on /i. From (35) we get 
(Va)^ = 5ji and from (36) we calculate {yiV^)j = (V2V3)* = 6ji {cj + idj) where 



Ki — K3 

Ki + K,3 



7 



- ft^3 



7 + iv, 



(a + /x)' + ^2' 



V — 



2/3^ 



(a + /i)' + /32' 



and 



«i (7^ - ^'^) + 267^^, 
2aj7z/ - 6 (7^ - . 



(A.5) 
(A.6) 

(A.7) 
(A.8) 



The general wobble tau functions are then given by 



r,- = 1 + 26" 



ttj cos C + & sin C - [a] + 6^) e"* 



+ 25,-ie''+'''' 



2a2 

/3 



+ 



9 cos C - sin C - ^ [a] + 6^) [c] + df) + -e"" 



(A.9) 
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Appendix A.l. Gardner hierarchy solutions 

The breather or wobble of the Gardner equation (25) is given by 

v = d^\iiy (A. 10) 

replacing (A.4) or (A. 9), respectively. Prom (45) we have 

77 = 2q;x + 2q; (a^ - 3^^) (A.ll) 
C =2px-2p{p^ -3a^)t, (A. 12) 

r]^ = 2fix + 2fj,H. (A. 13) 

For (27) the only modification comes from (46) that yields 

ri =2Q.x + 2a{a^ -10a^l3^ + 5/3^)t, (A.14) 
C ^2/3x + 2/3 {/3^ -lOa'^P^ + 5a^)t, (A. 15) 

rii,^2pLx + 2pLH. (A.16) 

Appendix A. 2. mKdV hierarchy solutions 

For the mKdV equation (17), the breather or wobble is given by 

v^vo + d^\n- (A.17) 

using (A.4) or (A. 9), respectively. In the coefficients (A.l), (A. 2), (A. 7) and (A. 8) we 
should replace /i — )■ wq- 

Taking into account the dispersion (49) we have 

77 = 2«x + 2a (a^ - 3/32 _ 1^2^ (A.18) 

C ^2(5x-2(5{(5'' -2,0" + lvl)t, (A. 19) 

?7„o = 2vqx - vlt. (A.20) 

For (18) the dispersion (50) imphes 

77 ^2ax + 2a {a^ - lOa^p^ + 5^^ - f u^a^ + fv^P^ + fv^) t, (A.21) 

C ^2Px + 2P (V - lOa^^^ + + f^^o'/S' - f^o«' + f ^0) (A.22) 

77^0 = 2voX + fv^i. (A.23) 
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